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In the effective-mass approximation, the step-like crystal potential of a wurtzite semiconductor heterostucture 
should be supplemented by Dirac delta-function heterointerface terms. They stem from the difference in the 
Bloch functions of the semiconductors and remain finite even for structures with graded chemical composition, 
where the terms are presented by a smeared Dirac delta function. We find these heterointerface potentials by 
employing the k ■ p method, and evaluate their strength from band-structure parameters of bulk materials. 
These potentials are weak for semiconductors compliant with the cubic approximation, which forces the zinc- 
blende crystal symmetry upon the wurtzite lattice. Nevertheless, they can produce a noteworthy effect due 
to a strong built-in electric field usually present in wurtzite heterostructures. We estimate that for GaN/AlN 
[0001] heterojunctions their net contributions to the energy of conduction and valence band states are 3 meV 
and 10 meV, respectively. The presence of the interface potential can modify the shape of the valence-band 
spectrum calculated without the potential. 
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I. INTRODUCTION 

In the last two decades, much effort has been de¬ 
voted to accurate determination of the band offsets 
and polarization fields in group 111-nitride semiconduc¬ 
tor heterostructures These material characteristics 
along with parameters of the bulk semiconductors are 
crucial elements of the band structure engineering of 
modern electronic and optoelectronic devices.—^— Mean¬ 
while, it is tacitly supposed that the knowledge of the 
band offsets and polarization fields provides rather accu¬ 
rate description of the effective potential of a heterojunc¬ 
tion. Gan the simple model of a step-like crystal potential 
be defective? 

Zhu and Kroemer were apparently the first to give a 
detailed analysis of the potential of an abrupt heterojunc¬ 
tion to be used in the effective-mass (EM) equations 
They showed that the heterojunction can be described 
by a potential step to represent the band offset and a 
heterointerface potential in the form of a Dirac delta 
function, the strength of which is subject to the details 
of the heterointerface on an atomic scale. This result 
of the tight-binding model has later been confirmed by 
consistent application of the k ■ p method )^^'^^ with the 
only input quantity, other than parameters of bulk ma¬ 
terials, being a phenomenological coordinate-dependent 
form-factor G(r), see Fig.[I] that characterizes the chem¬ 
ical composition of a heterostructure<^ 

Contrary to what one might expect, the interface po¬ 
tential of a hetero junction is not formally represented by 
a space derivative of the function G(r). Briefly, its origin 
can be traced as follows. The form-factor can be decom¬ 
posed into a gentle part, processed by using the standard 
technique of the EM approximation for slowly varying 
perturbations^!^ and a rapidly varying part, the space 
extent of which d is supposed to be small on a scale of the 



FIG. 1. Phenomenological one-dimensional form-factor G{z) 
of a single heterojunction. The transition region of the junc¬ 
tion is d, so that G(a)|^<_d /2 = 0 and G(«)|^><j /2 = 1- 


wavelength L of the envelope function (EE). The rapidly 
varying part can be processed in the spirit of the ap¬ 
proach presented by Bir and Pikua^^ (Chap. 27) for the 
central-cell corrections in the problem of shallow impu¬ 
rity levels in bulk semiconductors. This rapidly varying 
part can accordingly be formulated as a convergent se¬ 
ries of interface-localized generalized functions each 
successive term contributing less than the previous one in 
the parameter a/L, where a is the lattice constant. It is 
the leading term of the series that gives rise to the delta- 
function potential of Zhu and Kroemer. The strength of 
the interface potentials not only depend on the micro¬ 
scopic details of the heterointerface but also on the ori¬ 
entation of the structure’s growth direction with respect 
to the crystal axes.— 

Zhu and Kroemer pointed out that the heterointerface 
terms were small or negligible for GaAa/AlGaAs het¬ 
erojunctions but might be strong in other semiconductor 
systems . Since then, there have been obtained no ex¬ 
perimental evidences of the existence of the diagonal in 
the band index interface potentials, which couple states 
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of the same type. It is not surprising because even if 
they are present and properly included into the EM equa¬ 
tions, the potentials do not modify symmetry but merely 
spectrum of the system, which may have been misinter¬ 
preted in terms of adjusted EM parameters or band off¬ 
sets. On the other hand, the interface potentials can 
also produce interband couplings, which, though gener¬ 
ally weaker than the intraband ones, are easily identified 
by unique effects they produce. Thus it has been shown 
for [001]-grown zinc-blende heterostructures that the in¬ 
terface terms are responsible for such phenomena as the 
intervalley and X^-X^j^ mixing of the con¬ 

duction band states, mixing of light and heavy holes at 
the center of the two-dimensional Brillouin zone;^ and 
also anisotropy of the interband optical transitions^i^^ 
and a twofold anisotropy of the in-plane spectrum of 
the conduction band electrons in asymmetric quantum 
wells 

Here we will ignore the rapidly varying part of the 
function G'(r) to focus on another type of the interface 
potential, not considered so far, that is missing in zinc- 
blende heterostructures but present in wurtzite ones. It is 
proportional to a derivative of the form-factor G(r) and 
is thus finite even for structures with graded chemical 
composition where the rapidly varying part of the func¬ 
tion G(r) produces an exponentially small contribution^^ 
(Chap. 22). 

II. HETEROINTERFACE POTENTIALS IN THE k-p 
METHOD 


the same lattice constant. The strain-induced piezoelec¬ 
tric field and the spontaneous polarization field are taken 
into account via the inclusion of an external scalar po¬ 
tential W = W{t), while the strain modification of the 
band edge energies and the band offsets can be accounted 
for by using the strain Hamiltonian j^^i^^ For brevity, the 
latter as well as the spin-orbit interaction are considered 
only schematically. The strain-induced change in the ma¬ 
trix elements of the momentum operator p = (px,Py,Pz) 
between Bloch states, being small in the parameter of the 
order of strain (several percent typically), is neglected. 

Dealing with states near the Brillouin zone center, it 
is convenient to use the complete set of Kohn-Luttinger 
functions^ with the unit-cell orthonormalized periodic 
parts Uno = Uno (r) = 1''^); which are chosen real, specified 
as follows: 

T ~ (2) 

where toq is the free electron mass and e„o is the band 
edge energy of the nth band for the potential well mate¬ 
rial. We also introduce here zone-center Bloch functions 
Uno = Uno (r) = \n) and the band edge energies ?„o for 
the barrier material: 

T U2^ ‘^nO ~ CnO'^nO- (3) 

Treating 6U as a perturbation^ we have the following 
approximate expressions for the band edge energies: 

CnO — “f (4) 


Let US consider a heterostructure composed of two 
semiconductors, so that the single-particle potential en¬ 
ergy of an electron is 

U = U{r) = Ui+GSU, (1) 

where Ui = Ui (r) is the periodic lattice potential of the 
nominally basis (potential well) material, while SU = 
SU(r) = U 2 — U 1 is the perturbative periodic potential 
and U 2 = G 2 (r) is the periodic potential of the barrier 
material. We suppose that the phenomenological func¬ 
tion G = G (r) takes values of the order of unity or less 
even at the heterointerfacesIdeally, G can even be a 
step-like function so that G = 0 in the region of the po¬ 
tential well material and G = I in the region of the bar¬ 
rier material;^ We also suppose that the semiconductors 
composing the structure are related, which means that 
5U can be treated as a small perturbation as compared 
to the basis potential Ui. 

We do not consider in details deformation effects ap¬ 
pearing in lattice-mismatched pairs, which would lead 
to redundant complications, not essential for our results. 
We suppose that an inhomogeneous coordinate transfor¬ 
mation restoring the unstrained lattice of the basis ma¬ 
terial and adjusting the lattice of the barrier material to 
that of the basis one has been performed^ so that in 
the transformed space both materials are described by 


and for the band edge Bloch functions: 


‘^nO — '^nO T 


SUr,'r. 


n'^n 


^nO ~ ^n'O 




( 5 ) 


where SUn'n = {n'\dU\n). 

Let us suppose that Cartesian coordinates are oriented 
so that the coordinate 0 is along [0001] direction (c-axis) 
of wurtzite. Here we explicitly introduce the valence 
band Bloch functions Uxo and Uyo, which transform as 
the coordinates x and y, respectively, belonging to the 
representation Tg of the space group Ggu, and also the 
valence band Bloch function and the conduction band 
Bloch function Uco, which transform as the coordinate z, 
belonging to the representation Ti 

In the mean field approximation, the Schrodinger equa¬ 
tion readsi^ 

(r) = (r), (6) 

with the term Use comprising the spin-orbit interaction 
and the strain Hamiltonian.^?® Defining the nth band EE 
in r-representation as = Fn (r), we have the following 
expansion for the total wave function vk (r):^^ 


(r) = ^ 

n 


( 7 ) 
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where the summation is over all bands. The EFs satisfy 

I 


the following k ■ p system of equations in the coordinate 
representation: 





P ■ {n\p\n') 


Too 


G (r) SUn 


H. 


(se) 


Fn' = 0. 


( 8 ) 


Equation (|51) is obtained by treating the functions G (r) 
and W (r) as gentlei^i^i At this point, we neglect all 
central-cell like corrections due to a rapid variation of the 
function G (r) at the heterointerfaces .^^ As in Eq. m, 
the term (r) takes into account contri¬ 

butions from the spin-orbit interaction and the strain 
Hamiltonian, now acting on EFs. The system of Equa¬ 
tions ([8]) is valid for slowly varying functions 

The conduction and valence band EM equations can 
be obtained from the set of Eq. (jS)) by using the Lowdin 
perturbation schem o^^’^^ up to the second order. For the 
conduction band, neglecting the spin-orbit interaction, 
we have: 


(Efc + i?ce + AC/c)Fc = eFc, (9) 

where i?c is the ordinary conduction band Hamiltonian 
for a wurtzite heterostructure. It reads as follows: 

= Cco + + G5U,, + W, (10) 

zrriz ^rrix 

where the parameter 6Ucc defines the conduction band 
offset without strain contributions, ruz and rux are the 
conduction band EMs, 



with j = z,x, and iJce is the conduction band strain 
Hamiltonian) ^^:^^ 

Flee = 0,l£zz + 02 {Exx + £yy) , (12) 

where oi = ai(r) and 02 = 02 (r) are the material- 
dependent conduction band deformation potentials, and 
e = e(r) is the position-dependent strain tensor. 

The term AUc = A17c(r) in the left-hand side of Eq. ([H]) 
is the new heterointerface potential introduced in this 
work, 


AC/c(r) = S • (VG(r)), 
where the real vector S is 

g _ {c\p\n)SU„c 

imo ^ eco - Cno 

n^c 


(13) 


(14) 


It has been derived by Leibler^SS, but misinterpreted as 
producing a quasi-momentum shift in the position of 


the spectrum extrema, and thus concluded to vanish for 
high symmetry points of the Brillouin zone. Indeed, it 
takes place for T point states in cubic crystals but not 
in wurtzite ones. It follows from the symmetry argu¬ 
ments, due to the fact that the operator pz belongs to 
the representation T 1 of the space group Cgy that the 
Sz component of the vector S is still finite. 

For the valence band, described by the band’s EFs F^, 
Fy and F^, the system of the EM equations can be pre¬ 
sented as follows: 

E„/ = eF„, n = X, y, z. (15) 

n'=x,y,z 


The matrix is the ordinary valence band Hamil¬ 

tonian for wurtzite heterostructures, which includes the 
kinetic energy, step-like potential energy, spin-orbit in¬ 
teraction and strain Hamiltonians given in Ap¬ 

pendix!^ for reference. The matrix AU = AU(r) has 
the following finite components: 


AHxx = AHvv = Ex 


9G(r) 

^ dz ’ 


AUzz = K 


AGxz = -6NsG{v)px + Ex 


aG(r) 

’ dx ’ 


AC/,x = AUl = 5N^G{v)px + E 


dG{v) 


AUyz = -6N3G{r)py + E 


dG{r) 

' Ck ' 

dy 


AUzy = AUl = 6N3G{r)py + E: 


dG{r) 

dy 


dG{r) 
dz ’ 


(16) 


The binary material parameters entering the above ex¬ 
pressions are 


h {j\Pz\n)dUnj 

imo ^ EjO - EnO 


(17) 



h{x\px\n)dUnz 


(18) 


(19) 


2imo 
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E 

n^x,z 


1 


1 


^xO ^nO ^zO ^nO 

h{z\p^\n)dUn^ 

2imo 


( 20 ) 


The matrix AU has contributions of two types. The 
terms proportional to the parameter SN^ provide position 
dependence for the linear in the momentum operator el¬ 
ements entering the kinetic energy matrix of Eq. (IA 6 I) . 
which are proportional to N^, see Appendix El The rest 
contributions, proportional to first partial derivatives of 
the form-factor G'(r), are finite only at heterointerfaces. 


Exo = CyO = Czo + ^CR ~ EzOi and introduce the band 
gap energy Eg = Eco — Czo and the Kane matrix elements 
Pi = —ih{c\pz\z)/mo and P2 = —ih{c\px\x)/mo. With 
these notations and the adopted approximation we have: 

^ ,, PlSUcz ,, P2^C/cz 

Jz — Vzz— p , Vxz — _ , 

Eg Eg ( 21 ) 


To be consistent with the approximation of zero crys¬ 
tal splitting, we may use P 2 ~ Pi so that all interface 
potentials are governed by a unique constant a: 


III. ESTIMATES 


a = 


PlSUcz 

Eg 


( 22 ) 


To estimate the above binary parameters given by 
Eps. (|14p , (|lTp . (jl9p . and. (uni), we restrict ourselves with 
the Kane model , considering only conduction and 
three valence bands. However, it should be borne in 
mind that the estimation may lack accuracy for actual 
wide band gap materials such as GaN and AIN. We fur¬ 
ther neglect the energy Acr of the crystal splitting i^^i^^ 


The conduction band interface potential of Eq. (US is 
thus: 

A(7e(r) = , (23) 

while the valence band matrix AU is 


AU = 


/ 0 
0 


\SN3Gir)px + 


0 

0 


5NoG{v)py + 


-5NoG{v)px + a 

-6NoG{Y)py p a 
.dG(v) 


a- 


dz 



(24) 


Let us now demonstrate how knowing only bulk pa¬ 
rameters of the materials composing the structure it is 
possible to extract the non-diagonal band offset SUcz and 
evaluate the parameter a. We recall the following explicit 
expression for the parameter A 7 entering the kinetic en¬ 
ergy Hamiltonian of Eq. (IA 6 I) , see Appendix 


A 7 = 


h(zlpxlx} 


i\/2mo 

for the potential well material, and 

h(zlpx |x) 


A 7 — 


i\/2mo 


(25) 


(26) 


for the barrier material. Using Eq. dSD, adapted to the 
Kane model, we have: 


From Eqs. 


(z| ~ (z| - |x) ~ |x). (27) 


and ( 1221 ) we immediately obtain: 


6U,z = 


V2Eg (At - At) 


P 2 


and, considering P 2 ~ Pi, we finally have: 




At — Ai 


(28) 

(29) 


The strength of the interface potentials is very weak 
for systems based on wurtzite materials compliant with 
the cubic approximation's Eor such semiconductors, in 
particular |Pi| Ri IP 2 I |At|. For example, using the 
parameters of GaN^S mz = 0.2mo and Eg = 3.5 eV, we 
have Pi ~ 8 eVA, while At ~ 94 meVAiSiSS Vurgaft- 
man and Meyer— report At = 0 for wurtzite AIN, so that 
for GaN/AlN heterostructures a r; 130 meVA, while the 
interband offset 6Ucz is as small as 60 meV. The latter 
estimate has already been reported— S Under these con¬ 
ditions, the inclusion of the interface potentials in the 
EM Hamiltonians cannot modify the states essentially. 
In particular, one may not expect that they can produce 
Tamm-like states localized at heterointerfaces. 


Let us evaluate the effect of the interface potential on 
the spectrum of the conduction band states in a quantum 
well formed by a single heterojunction GaN/AlN in the 
model of the mathematically abrupt form-factor G(z) = 
Q{z), where 0(z) is the Heaviside step function, so that 
d<d{z)/dz = (5(z), see Fig.O where 6{z) is the Dirac delta 
function. Neglecting for brevity the strain Hamiltonian 
Pee in Eq. we consider the Hamiltonian = 
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FIG. 2. Potential energy profile of the conduction band for 
a heterojunction GaN/AlN, solid line, including the interface 
potential, arrowed line. Also shown are the ground eigenen- 
ergy, dotted line, and the corresponding envelope function, 
dashed line. 

Hc + AUc-. 

= 6^0 + :^ + + 0(z),5t/cc + W{z) + a5{z). 

(30) 

As soon as the term AUc = ad{z) is a small perturba¬ 
tion, the discrete spectrum of the Hamiltonian of Eq. o 
can be approximated as follows: 

= +a{N\S{z)\N), (31) 

where are eigenenergies of the Hamiltonian He with 
the corresponding EFs = |Af), that is Hc\N) = 

Diagonal matrix elements of the commutator [pz,He\ 
disappear for any two-dimensional subband index N: 

{N\pzHe-HePz\N)=0, (32) 

resulting in the following identity)^ 

SUec{N\S{z)\N) = -(iV|^ffi|iV). (33) 

dz 

The right-hand side of Eq. (1551) is proportional to the 
average build-in electric field mostly originating due to 
the spontaneous polarization.^? If we neglect the con¬ 
tribution of the AIN barrier region to the matrix ele¬ 
ment on the right-hand side of Eq. (1551) . we will have: 

Ri —eifcaN, where —e < 0 is the electron 
charge and EcaN is the polarization field in the quan¬ 
tum well material (GaN). If we use Slice = 2.5 eV, which 
is accepted as the conduction band offset for a hetero¬ 
junction formed by AIN lattice matched to GaN^i^ and 
EcaN = 5 MV/cm,— we have the following estimate: 
(iV|5(z)|Af) « 0.02 A-i, and a{N\S{z)\N) « 3 meV for 
the values of the parameters Ar given by Vurgaftman and 
Meyer 

For the valence band, the above arguments expressed 
by Eqs. dSlD and (1551) are not directly applicable be¬ 
cause the interface potential is not purely diagonal. As 


usual, non-diagonal components of the potential make a 
minor contribution, but the potential still has a diago¬ 
nal component to influence the ^-component of the EF. 
Gonsequently, when the interface potential is included, 
the spectrum of the hole states may not be uniformly 
shifted as we observed in the conduction-band case, but 
the in-plane effective masses of the subbands will also be 
modified. The interface potential has the strongest influ¬ 
ence on the states with nearly isotropic in-plane electron 
density. For such states, the conduction-band estimation 
approach can be used. Then, for a heterojunction formed 
by GaN lattice matched to AIN, we must take into ac¬ 
count smaller band offset SUzz = 0.85 eV— which results 
in the contribution of the interface potentials of the order 
of 10 meV. 

Note that besides the restrictions of the four-band 
Kane model, applicability of which to wide band gap 
semiconductors is questionable, the values of the bulk 
parameters A 7 cannot be considered as established in 
the literature, so that the above figure for a may thus 
be unreliable. For example, using the parameters given 
by Rinke et al.^ Ay m 46 meVA (GaN) and A 7 r: 
26 meVA (AIN), we derive a R 30 meVA, which is four 
times as small as the figure based on the parameters given 
by Vurgaftman and Meyer 


IV. CONCLUSIONS 

Using the k ■ p method, we have shown that the po¬ 
tential of a wurtzite semiconductor heterostructure must 
be supplemented by Dirac delta-function terms located 
at heterointerfaces. These heterointerface potentials are 
present even for structures with graded chemical com¬ 
position, and their strength is governed by binary ma¬ 
terial parameters that can be deduced knowing only 
band structure characteristics of the bulk semiconduc¬ 
tors. They are weak for materials whose wurtzite lattice 
slightly deviates from that of a cubic crystal. We have 
obtained these potentials for electrons and holes and es¬ 
timated their strength. In the four-band Kane model, 
the potentials are defined by a unique parameter pro¬ 
portional to the difference in the parameters A 7 of the 
semiconductors of the structure. For GaN/AlN [0001] 
heterojunctions, we have evaluated that the contribution 
of the potentials to the energy of conduction and valence 
band states is 3 meV and 10 meV, respectively. In addi¬ 
tion, the presence of the interface potential can modify 
the shape of the valence-band spectrum. 
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Appendix A: Valence band Hamiltonian 


In the basis of the band edge Bloch functions Ux, Uy 
and Uz, the valence band Hamiltonian in the EM approx¬ 


imation without interface potentials is as follows: 

H(v) (A1) 

where is the position-dependent potential energy of 
an electron: 


H(o) 


/cxo + G<5I/xx + W 0 0 

0 CyO + GSUyy +W 0 

y 0 0 Czo + GSUzz + W 


(A 2 ) 


where Cxo = CyO) while SU^x = ^Uyy and SU^z are the valence band offsets, not including stain effects, which appear 
in the strain Hamiltonian 


h(=) 


( ^l^xx -\- Tll2^zz ’^l^xy 

'^l^xy '^l^xx “t ^l^yy “f ‘^2^zz 

^2^xz ^2^yz 


‘^2^xz 

‘^2^yz 

^3 (y^xx “ 1 “ ^yy ) + h£ ZZ 


(A3) 


where real material-dependent parameters /i, I 2 , mi, m 2 , 
ni, and n 2 are expressed via conventiona l^'^^i^^'^^ com¬ 
ponents of the deformation potential tensor as follows: 

h = D2 -l- D4 + Z? 5 , mi = D2 + D4 — D5, 

ni = 2D5, I2 = Di, m2 = Di D3, (A 4 ) 

n2 = V2 Dq, m 3 = £> 2 - 


where ax, ay, and CTz are the Pauli matriceSf^l and real 
A 2 = A 2 (r) and A 3 = A 3 (r) are the parameters of the 
valence-band spin-orbit splittingj ^^i^^ The valence band 
EFs Fx, Fy, and Fz are thus defined as two-component 
spinors. 


The Hamiltonian in Eq. IMI describes the spin- 
orbit interaction i^^i^^ 

( 0 -i/S.2az iAsay \ 

iA 2 az 0 —iA^ax ] , (AS) The kinetic energy Hamiltonian of Eq. (lAll) is 

-iA^ay iA^ax 0 /| 


/ LipI + MipI + M 2 PI NiPxPy N2 PxPz - Nspx \ 

= NipxPy MipI + LipI + N 2 PyPz - Nspy , (A6) 

\ N2PxPz + NsPx N2PyPz + NsPy Ms {pi + p^ -H £2Pz/ 


where 
Li = 


{A 2 -I- A 4 -I- As), L 2 = - —Ai, 

2mo 


2 mo 

Ml = -— (A 2 -I- A 4 — As), Ni = f, 2As, 


2 mo 


2 mo 


(A7) 


M2 = -— (Ai -I- A3), N2 = - —V^Ao 
Zttiq Zttlq 

Ms = -—A 2 , Ns = -V^Ar, 

2mo n 


(AS) 


and Ai, A 2 , ... A 7 are real material parameters in con¬ 
ventional notations.— 

If it is preferable to have the valence band EM equa¬ 
tions expressed in terms of the EFs Fi, F 2 , and F 3 


that multiply the conventiona P^i^^ basis functions ui = 
{ux + iUy)j'/2, U2 = (ux — iuy)l\/2, and us = iuz, one 
can use the canonical transformation 



so that Eq. m now reads 

^ = eF„, n = 1, 2, 3, 

(AlO) 

where = T-^H^T and Atj == T-^AUT. 
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